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It is shown that the quantum massive non-Abelian field theory established in the former papers is 
renormalizable. This conclusion is achieved with the aid of the Ward-Takahashi identities satisfied by 
the generating functionals which were derived in the preceding paper based on the BRST-symmetry 
of the theory. By the use of the Ward-Takahashi identity, it is proved that the divergences occurring 
in the perturbative calculations for the massive gauge field theory can be eliminated by introducing 
a finite number of counterterms in the effective action. As a result of the proof, it is found that 
the renormalization constants for the massive gauge field theory comply with the same Slavnov- 
Taylor identity as that for the massless gauge field theory. The latter identity is re-derived from 
the Ward-Takahashi identities satisfied by the gluon proper vertices and their renormalization. To 
illustrate the renormahzabihty of the theory, the renormalization of the QCD with massive gluons is 
, . concretely performed and the one-loop effective coupling constant of the QCD is derived and given 

' an exact expression in the renormalization. 
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lO . I- INTRODUCTION 

> _ 

■ As was mentioned in our previous paper^^l , in the original attempts of setting up the massive non-Abelian gauge field 
theory without Higgs mechanisml^"^"] , the massive Yang-Mills Lagrangian itself was considered to form a complete 
formulation of the massive gauge field dynamics and used to establish the quantum theory. In the theory, there are 
two problems which were announced to be difficult to solve: one is the gauge-non-invariance of the mass term in 

OO , the action, another is the nonrenormalizability of the quantum theory. In Ref. [1], the first problem has readily 
been circumvented from the viewpoint that the massive gauge field only exists in the physical space spanned by the 
' transverse part of the vector potential. In this space, the action of the massive non-Abelian gauge fields whose masses 

■ are taken to be the same is gauge-invariant with respect to the infinitesimal gauge transformations which are only 
necessary to be taken into account in the physical space . If we want to represent the massive gauge field dynamics in 

D the whole space of the vector potential, the massive gauge field must be viewed as a constrained system. The Lorentz 
' gauge condition, acting as a constraint, must be introduced initially and imposed on the Lagrangian expressed by 
J> . the full vector potential. From this point of view, it has been shown that the massive gauge field theory can well be 
established on the basis of gauge invariance. In the preceding paper (referred to as paper I later on). It was shown 
. , that the quantum massive gauge field theory described in Ref. [1] has a BRST-symmetry, that is to say, the effective 

■ action and the generating functional of Green functions are invariant with respect to a kind of BRST-transformations. 
" ■ ' From the BRST-symmetry of the theory, a series of Ward-Takahashi (W-T) identitiesl^""^^' obeyed by the generating 

functionals were derived and used to prove the unitarity of the theory. 

In this paper, we are devoted to proving the renormahzabihty of the quantum massive gauge field theory. The 
renormahzabihty of such a theory may be seen from the intuitive observation that the Feynman rules derived from 
the effective action, except for the gluon and ghost particle propagators, are the same as those given in the massless 
gauge field theory, and the massive propagators have the same ultraviolet behavior as the massless ones. In particular, 
the primitively divergent diagrams are completely the same in the both theories. This fact suggests that the power 
counting argument of demonstrating the renormahzabihty is useful not only for the massless gauge field theory, but 
also for the massive gauge field theory. Theoretically, to accomplish a rigorous proof of the renormahzabihty of 
the massive non-Abelian gauge field theory, it is necessary to implement a subtraction procedure to see whether 
the divergences occurring in perturbative calculations of Green functions and S-matrix elements can be removed by 
introducing a finite number of counterterms in the action . This procedure, as one knows, amounts to the well-known 
R-opcration invented by Bogoliubov, Parasiuk, Hepp and Zimmermann^^'^"^^!. The principal idea of proving the 
renormahzabihty of the theory under consideration is the usage of the W-T identities derived in paper I. In the 
Landau gauge, as mentioned before, these identities formally are identical to those for the massless gauge field theory. 
Therefore, the proof of the renormahzabihty almost is the same as for the massless theory. From the proof, it will 
be seen that the divergences occurring in perturbative calculations can surely be eliminated by introducing a finite 
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number of counterterms. As a consequence of the proof, it will be found that the Slavnov- Taylor identityl^^^^^'l for 
the rcnormalization constants which was derived in the massless gauge field theory also holds for the massive gauge 
field theory. As will be shown, the Slavnov- Taylor identity can also be derived from the W-T identities satisfied 
by the vertices and their renormalization. To illustrate concretely the renormaliz ability of the theory, the one-loop 
renormalization of the quantum chromodynamicm with massive gluons (will be called massive QCD) will be performed 
and the one-loop effective coupling constant for the QCD will be derived and given an exact expression. From this 
expression, the previous result for the QCD with massless gluons which was obtained by the minimal subtraction'^** 
will be recovered in the large momentum limit. 

The remainder of this paper is arranged as follows. In section 2, we describe the general proof of the renormalizability 
of the massive non-Abelian gauge field theory. In the proof, we start from the generating functional of proper vertices 
and use the W-T identity for this generating functional to show that the divergences occurring in the perturbative 
expansion of the generating functional can completely be eliminated by introducing a finite number of counterterms 
in the action. The counterterms are explicitly given by solving the equations obeyed by the counterterms which 
are derived from the W-T identities. In section 3, we analyze the role played by the counterterms and prove the 
equivalence between the additional renormalization and the multiplicative renormalization. as a result of the proof, 
the Slavnov- Taylor identity satisfied by the renormalization constants is naturally derived. Section 4 is used to derive 
the W-T identity satisfied by the gluon three-line proper vertex, discuss its renormalization and show how the identity 
for the renormalization constants can follow from the above-mentioned W-T identity. In section 5, the same things 
will be done for the gluon four-line proper vertex. The last section serves to make some concluding remarks. Section 
6 serves to derive the one-loop effective coupling constant for the QCD with massive gluons. In Appendix, we quote 
the proof which shows how the counterterms written in section 2 satisfy the W-T identities. 



II. W-T IDENTITY AND COUNTERTERMS 



From this section, we start to prove the renormalizability of the massive non-Abelian field theory in which all 
the fields have the same masses. The theory concerned typifies the quantum chromo dynamics (QCD) with massive 
gluons'^^l, the SU(2)-symmetric hadrodynamics^ or some others. In the following, we will take the QCD with massive 
gluons as an example. Since all of vertices and even Green functions in the theory may be derived from the generating 
functional T for proper vertices which was defined in Eq. (3.11) in paper I, it is only necessary to deal with the 
renormalization of such a generating functional and to see whether the divergences appearing in the perturbative 
expansion of the functional F can be removed by introducing a finite number of counterterms . 

Let us make use of the loop diagram expansion, which is a power series in the Planck constant h, for the proper 
vertex generating functionall^*^'^^! 

oo 

f = ^f„ (2.1) 

where F is, according to Eq. (3.21) in paper I, defined by excluding the gauge-fixing term 

t = T+^J d^d^Alf. (2.2) 

In the tree diagram approximation, as one knows, the proper vertex generating functional Fq just is the generalized 
action which is defined by the ordinary effective action plus the BRST-external source terms as was shown in Eq. 
(3.4) in paper I and, according to the definition given in Eq. (2.2), it will be rewritten as 

To = So=to-^J d^xid^A-^y (2.3) 

where 

fo = 5*0 = / d^x{i;[i^'^{df, - i.g^;i7^")]V' - M^^P - {F'^^'^F^, 

+\w?A''^Al + C''d^{Vf&) + (Alp + AtPC + m^'^AA;^ (2.4) 
+v"-AC''}. 

In the above, represent the gauge fields, ip and ip denote the quark fields, m and M are the gluon and quark 
masses respectively. 
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Vf{x) = !^dl + Df{x) (2.5) 

here jj^ = arn? and 

Df = 6''% - gr'^Al (2.6) 

which is the ordinary covariant derivative, A'^, A^^, Aj4^ and AC" are the composite field functions which are defined 
from the BRST-transformations by dropping out the infinitesimal Grassmann number ^ and u"'', u", C, and C are the 
corresponding BRST-sources. 

It should be emphasized that according to the additional renormalization scheme which is adopted in the latter 

proof, all the field fimctions and parameters in the effective action shown in Eqs. (2.3) - (2.6) arc rcnormalized 
ones and hence finite. In pcrtm'bativc calculations, each loop term in Eq. (2.1) is divergent and therefore has to be 
regularized by an appropriate regularization scheme which must be chosen to preserve the BRST-symmetry of the 
theory. Suppose the n-th term r„ in Eq. (2.1) has been separated into a finite part and a divergent part 
through the regularization procedure 

f„-r^ + f^ (2.7) 

In perturbative calculations, the divergences included in Eq. (2.1) may be eliminated order by order through a 
recursive construction of counterterms in the action. For instance, to eliminate the one-loop divergence Ff which is 
generated by using the action shown in Eq. (2.4) in the following perturbative expansion 

f[^o] =fo[^o]+fi[^o] + = 5o+ff[So]+f([5o] + , (2.8) 

we may choose a counterterm AS'o such that 

A5o = -ti[So] (2.9) 
whose concrete form will be given later. It is apparent that when we use the action 

5i = So + A5o (2.10) 

to recalculate the functional V of order Ti 

f[5i] = f[5o + A5o] = fo[5o + A5o] + fi[5o + A5o] + 

= 5o + A5o + ff[So] + ff[5o] + (2.11) 

= 5o + ff[5o] + , 

the divergence Ff [^o] in it disappears. In general, to remove the divergent part F^ in the n-loop term r„ of order 
we need to introduce a counterterm like this 

A5„_i = -ti. (2.12) 

Adding it to the action, we have 

Sn = Sn-1 + ASn-l (2.13) 

where the Sn-i has included the counterterms up to the order h^~^. The action Sn used to calculate the r„ will lead 
to a finite result in the n-th order perturbation. 

The counterterms mentioned above may be determined with the help of the following W-T identity satisfied by the 
generating functional of proper vertices which was derived in paper I 

f*f + a; = (2.14) 

and the ghost equation 



3 



where 



^ p_ f 4^ St ST 6f ST ST6T 6T ST 

J ^^SA^^S^^SC^S^^Si^Jc^S^'SC^ ^ ^ 

and 

w = y d^xd^AlixJCix) (2.17) 

As mentioned in paper I, in the Landau gauge, since /j. = and iv = due to i9'^A^(a;) = 0, the above W-T identity 
and ghost equation are identical to those for the massless gauge theory. In the loop expansion, when Eqs. (2.1) and 
(2.7) arc substituted in Eq. (2.14), one may obtain a scries of identities satisfied by the finite and divergent parts of 
the r„. Each of the identities contains terms which are of the same order of divergence. The identities of order h"' 



t(*Ti_,+ti_,*t(^o, 



(2.18) 



where 



plS) = — -— + — -— (2 19) 

Sipi Sui 6ui Sipi 

here (pi and Ui{i = A, ^l), tf), C) stand for the field variables A^, xj), ip, C" and source variables u^, C, (, u", respectively, 
and the symbol "." in each term on the right hand side (RHS) of Eq. (2.19) is an abbreviation notation of the 
integration as shown in Eq. (2.16). In the above, the fact that the uj is finite and of zeroth order of h has been 
noticed. This fact is obvious because we start from the renormalized field functions and parameters in the additional 
renormalization. 

Furthermore, the action constructed in Eq. (2.13) is also required to fulfill the W-T identity 

Sn* Sn+LU ^0. (2.20) 

On substituting Eq .(2.13) into Eq. (2.20) and noticing Eq. (2.12), one may find three equations[^°~^^l 

p(5„_i)f ^ = 0, (2.21) 

f n * = (2.22) 

and 

Sn-l * Sn-1 + w = (2.23) 

In Eq. (2.21), the operator p{Sn-i) is defined as 

p(S„-A - '^'^""^ ^ I '^'^""^ ^ (2.24) 
" Sipi Sui 6ui Sipi ' 

When we set the subscript n = 2, 3, • • • and repeatedly apply Eqs. (2.12) and (2.13), it is easy to find that the second 
term on the left hand side (LHS) of the last equation in Eq. (2.18) equals to zero. Thus, the last equation in Eq. 
(2.18) leads tol^^'"! 

p{So)ft = (2.25) 

where p(Sq) was defined in Eq. (2.19). In addition, when Eqs. (2.13) and (2.12) arc substituted into Eq. (2.15), one 
may obtain the ghost equations satisfied by the finite part TI_ and the divergent part FjJ of the functional r„. For 
the divergent F^, we have'^'^'^^l 
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(2.26) 



It is emphasized here that even in the general gauges {a ^ 0), the W-T identity and the ghost equation satisfied by 
the divergent functional Tf^ are the same as those given in the Landau gauge {a = 0) because the functions oj in Eq. 
(2.14) and the fi'^Cix) in Eq. (2.15) are finite. 

Clearly, the divergent functional FJ^ and thus the counterterm ASn-i may be determined by solving Eq. (2.25) 
and (2.26) together, or, instead, by solving Eq. (2.21) and (2.26) provided that the action S'„_i has been given in the 
former n — 1 steps of recursion. The general solution to the above equations was already found in the literaturel^'^'^^1. 
It consists of two parts as shown in the following 

ASn-l = -ti = J2 <Ha + p{Sn-l)Fn (2.27) 

a 

where the Ha in the first term are functionals of the field variables and tp which are invariant with respect to 

the gauge transformation. Therefore, they obviously satisfy Eq. (2.21) or (2.25). These functionals Ha can only come 
from the first four terms in Eq. (2.4) and of the forms 

^G = -\J d^x{d^Al - d^Al + gP'^^AlAlf, (2.28) 



Hf = j d^x^irid^ - igA;T-)^p, (2.29) 
Hm = j d'^xM^)^ (2.30) 



and 



Hm = J d^x{A^^^f (2.31) 

here Hm is written only for the transverse part of the vector potentials because only for this part of the vector 
potential, the action of mass term is gauge- invariant. The second term in Eq. (2.27) directly follows from the 
nilpotency property of the operator p, p^F^ = 0, with the functional F„ being arbitrary!^"! , as shown in Appendix. 

It is noted that the term ASn-i as a part of the action demands that the Fn must be a functional with minus mass 
dimension and minus ghost number so as to make the action to be dimensionless and of zero- ghost number . This 
can be seen from the expression 

piS^.,)F^ = J d x{^ • ^ + ^ • (2-32) 

by noticing that the ghost numbers of the functions A" , ip, ip, C°, C°, u^, u", C, C are 0, 0, 0, —1, +1, —1, —2, —1,-1 and 
the mass dimensions of these functions are 1, |, |, 1, 1, 2, 2, |, |, respectively. Furthermore, the F„, as easily seen, 
must satisfy the ghost equation. 

^_a„^ = 0. (2.33) 

With the requirements stated above, the form of the functional F„ will be uniquely determined, as given in the 
following 

F„ = / d^x{blA''^{ul - a^C") + J2 b':^^u,}. (2.34) 

It is noted that the coefficients in Eq. (2.27) and 6" in Eq. (2.34) all depend on the regularization parameter, say, 
the £ = 2 — n/2 (which tends to zero, when n ^ 4) in the dimensional rcgularizationl^^l. The operator p{Sn-i) in 
Eq. (2.32) implies that we have chosen the counterterm AS'^-i to be the solution of Eq. (2.21) for convenience of 
later recursion. 
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III. EQUIVALENCE BETWEEN THE ADDITIONAL AND MULTIPLICATIVE RENORMALIZATIONS 



Up to the present, the counterterm ASn-i appearing in Eq. (2.13) has exphcitly been constructed as given in Eqs. 
(2.27)-(2.31) and (2.34). The action Sn-i constructed in the foregoing steps of recursion has the same functional 
structure as that for the 5*0 given in Eq. (2.4). This can be seen from the fact that Eq. (2.21) has the same form as 
Eq. (2.25). An interesting thing is that the role played by the counterterm in Eq. (2.13) which has been found in Eq. 
(2.27) is only to make a change to the variables of the first term Sn-i in Eq. (2.13). If the coefficients in Eq. (2.34) 
are assumed to be infinitesimal, we have the following variations: 

OUi 

5F 

Sul-d,SC- = --^^,i = A. (3.1) 
According to the definition given in Eq. (2.24) and applying the ghost equation obeyed by the action Sn-i 

and the variations defined in Eq. (3.1), one can derive'^^"^^! 

Sn-l[(Pi,Ui] + p(Sn-l)Fn[ipi,Ui] 
— Q 1 _1_ '^■^"-1 SF„ I A\S„_i SF„ 

' \ sipi ' Sui Sui ' S(fii ' (3-3) 



i=A,il),il>,C 

= S„-i[iPi,u'i\ 



in which 



ip'i = (Pi+ Sfi = Y^(fi, i = A, ip^ip, C, 

v!^ = u, + hL, = Yf\u i = A,i;,ij,C, (3.4) 

where the coefficients Y^, according to the definitions in Eq. (3.1), can be calculated from the expression written in 
Eq. (2.34). The results are 

Y,^ = l + h\',i = A,i^,^l^,C. (3.5) 

Considering that the functionals in the first term of Eq. (2.27) are gauge-invariant and of the same functional structure 
as those terms in the Sn-i which are the functionals of the fields A^, ■)/) and ip, we are allowed to change the field 
variables from to ip^ in the functionals in the first term of Eq. (2.27) and combine these functionals with the 
corresponding terms in the S'„-i wj] shown in Eq. (3.3) together. The results is just to redefine the variables and 
physical constants in the action Sn-i- Thus, the action in Eq. (2.13) can be written as 

SnWi,Ui] = Sn-l[^/Zfvi, ^f^Ui] (3.6) 

where Z" and are the n-th order multiplicative renormalization constants for the field and source variables 
respectively. Eq. (3.6) establishes a recursive relation of the renormalization. When the order n tends to infinity, we 
obtain from Eq. (3.6) by recursion the following result 

SyuUi]=Soyi,u^] (3.7) 
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where 



'■p'i — \fZ'i'-pi; '"i' — vZiUi, 



g° = Zgg, = ZmM, m° = Z^m ^^"^^ 

are the bare quantities appearing in the unrenormahzed action Sq. The renormalization constants in Eq. (3.8) are 
defined by 

oo „ oo _ oo 

z,= u zr, z,= u zr, Zg=u z^, 

n—1 n—1 n—1 /o q\ 

oo oo W / 

Zm = n Zm, Zm = Yl Z^. 
n=l n=l 



Eq. (3.7) shows us that the ronormahzcd action has the same functional structure as the unrenormahzed one. In the 
derivation shown above, we start from the additional renormalization by introducing the counterterms to remove the 
divergences appearing in the functional T and then obtain the result of multiplicative renormalization which is given 
by introducing the renormalization constants to absorb the divergences and redefining the field functions and physical 
parameters. This just shows the equivalence between the both renormalizations. 

To be more specific, let us describe the one-loop renormalization of the functional Fi starting from the action 
written in Eqs. (2.3) and (2.4). For convenience of statement, we at first discuss the renormalization in the Landau 
gauge in which the gluon fields are transverse. For the one-loop renormalization, as indicated in Eqs. (2.9) and 
(2.10), we have to introduce a counterterm whose general form was given in Eq. (2.27) with the order label n = 1. 
In the first term of Eq. (2.27), the gauge-invariant functionals were written in Eqs .(2.28)-(2.31). The corresponding 
coefficients in Eq.(2.27) will be written as ao, ap, dM and Um- In the following, the order label will be suppressed and 
the source terms in Eq. (2.4) will be omitted for simplicity because these terms act only in the intermediate stages of 
the proof. As demonstrated in Eq. (3.3), the variables of the action Sq in Eq. (2.10) which was explicitly written in 
Eq. (2.4) will be changed to the ones shown in Eq.(3.4) owing to the effect of the counterterm given by the second 
term in Eq. (2.27) with the functional F being represented in Eq. (2.34) and the variables in the counterterms of the 
first term in Eq. (2.27) which are explicitly written in Eqs. (2.28)-(2.31) may also be made such changes due to the 
gauge-invariance of the functionals. Thus, the action in Eq. (2.10) without the source terms may be written as 

/ d^xlYFY^Y^^iji^id^ - igYAAlT-)^ - MYMYjY^^i^ , . 

-{YlYcid^Al - d,A- + gYAr^'^A^^Alf + \Y^Ylrn? A-^^ A^ ^ • ' 

+YX^YcC''UC°' + Ycgf^^di^CC^Al] 

where 



Yp = l + aF,YM = 1 + ttM, Yg = 1 + acYm = 1 + am (3.11) 

and the subscript T marking the transverse gauge fields A^ has been suppressed for simplicity of representation. 
When we define renormalization constants by[^^~^^l 



Z2 = YpY-^Y^, Z3 = YgYI = YcYx\ 

Z( = YpY^Y^^Ya, Z, = YgYI Z4 = YgYI (3.12) 

G 



Zi — Yc, Zm — YmYp ^, Z"^ — Y^Y^Z^ ^ — Y^Yq ^ 



and noticing the relation given in Eq.(2.3), We may write the full action as follows 

S = J d^xiZ^'^iijf'd^ - MZm)^ + gZ[i}rA''^T-^P 
-IZsid.A- - d.A^f - ^gZ.r'-id.A- - d,.A-)A''^A-'^ 

-\g'^Zif''^''f'"^''A^''A''''A'lAl + IZl^Z^w? A^A^t" ^ > 

-h.^aC'^DC" + gZif'''>'=di^C"-C''A''^ - ^(a''A°)2}. 

It is noted here that the last equality in Eq. (3.12) gives the definition of gluon mass renormalization constant as 

]^ /2 1/2 

Zm = Z^ ' Ym Ya which is consistent with that given in the renormalization of gluon propagator as shown in Eq. 
(4.25) of paper I. The action in Eq. (3.13) with the renormalization constants being defined in Eq. (3.12) just gives 
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the recursive relation shown in Eq. (3.6) with the label n = 1. This action would eliminate the divergence appearing 
in the generating functional T evaluated in the one-loop approximation. 
From Eq. (3.12), it is clear to see that 
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Zi Z\ z, 

Zi ^3 Z\ 



(3.14) 



This is the Slavnov- Taylor (S-T) identity satisfied by the renormalization constants which is completely the same as 
given in the masslcss QCDI^^''^^!. 

If we define the bare quantities as^^^'^^' 

Z3C^ CS = ^zC\ go = Z^Z-'/^g, (3.15) 



Mo = ZmM, mo = Zmiri, ao = Z^a 
and use the identity in Eq. (3.14), wo arrive at 

SiA;,^,^,^^^^] = 5o[Ar,V5o,^o,C'o",Co"] (3.16) 

where 

So = J d^x{Mh^.{d^ - igoATT-) - Mo]^o 

_ 1 (a^Ag- - a^A^^ + 5o/"''^A^M8'')2 , 
+CSDC§ + gor'^'d^CSC^Al^' ^''■"> 
+im2ArAg^-^(a,A»'')n 

is the unrenormalized action. Eqs. (2.3), (2.4), (3.16) and (3.17) indicate that the actions, renormalized and un- 
renormalized, have the same structure and thus the same symmetry, just as we met in the massless gauge field 
theory. 

We note here that although the above results are obtained in the one-loop renormalization. they c;an, actually, be 
considered to be the exact ones. In fact, for removing the two-loop divergence in the F, obviously, the second cycle of 
recursion of the renormalization can be carried out in the same way as stated in Eqs. (3.10)-(3.13) by starting from 
the action shown in Eq. (3.10) and higher order recursions can be continued further along the same line. The results 
given in Eqs. (3.10)-(3.17) formally remain unchanged for each cycle of the recursion. Therefore, by the recursive 
procedure, all the results denoted in Eqs.(3.12)-(3.17) can be regarded as the ones as shown in Eqs. (3.7)-(3.9). 

Now, let us describe the renormalizability of the theory in the general gauges. We firstly note that the results in the 
Landau gauge as given before can readily be extended to the other gauges. As mentioned before, in the general gauges, 
we still have the equations written in Eqs. (2.21) and (2.26) for the divergent part F^. Therefore, the counterterm, 
as the solution to the equations (2.21) and (2.26), is still expressed by Eqs. (2.27)-(2.31) and (2.34) with a note that 
except for the gluon mass counterterm written in Eq. (2.31) which is still given by the transverse field, the vector 
potential in the other counterterms now becomes the full one. It is clear that the statements in Eqs. (3.1)-(3.9) 
completely hold for the general gauges. The results described in Eqs. (3.10)-(3.17), as easily seen, except for a few 
supplements for the gluon and ghost particle mass terms, are also preserved in the present case. The gluon mass 
term in Eq. (2.4) is now written for the full vector potential. When the countertems in Eq. (2.27) with their explicit 
expressions denoted in Eqs. (2.28)-(2.31) and (2.34) are added to the action, the mass term becomes 

I d^x]Y^Ylm^\A%^A^^ + y-^A^Al^]. (3.18) 

This term should replace the corresponding term in Eq. (3.10) to appear in the action. In the above, the factor Y\ 
arises from the variable change generated by the counterterm given in the second term in Eq. (2.27) as described in 
Eqs .(3.3) and (3.4) and the factor comes from the counterterm denoted in Eq. (2.31) which is written only for 
the transverse fields . When we notice the last equality in Eq. (3.12) and define 

^3 = (3-19) 

and 

A"^ = A^" + V^^r- (3-20) 
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Eq. (3.18) can be written as 

\ I d'xZiZsm'{A^/A^^^ + Z'^A^Al^) = d^xZ^ZWA'^^A; (3.21) 

where the orthonormahty between the transverse and longitudinal vector potentials has been considered. With the 
above expression, the gluon mass term in Eq. (3.13) can be understood for the full vector potential in the general gauge. 
Eq. (3.21) shows that the renormalization of the longitudinal part of the vector potential is different from that for the 
transverse part by an extra renormalization constant \/Z^- This result is consistent with the renormalization of gluon 
propagator shown in section 4 of paper 1 where the longitudinal part of the propagator has an extra renormalization 
constant Z'^. According to the definition given in Eq. (3.15), the gluon mass term in Eq. (3.21) will come to the form 
expressed by the bare quantities as shown in Eq. (3.17). 

Let us turn to the renormalization of the ghost particle mass term. As mentioned in section 2, the ghost particle 
mass term /i^C" in Eq. (2.15) which appears in general gauges is finite and therefore it is not included in the ghost 
equation (2.26) satisfied by the counterterm F^. That is to say, the ghost particle term C°'ji'^C'^ in the Lagrangian 
remains unchanged in the process of divergence subtraction because it is not affected by the counterterm. The 
counterterm only changes the kinetic energy term C"'nC"' and the interaction term gf'^^'^d^C'^C^A'jj^ as for the case 
in the Landau gauge. This fact coincides with the renormalization of ghost particle propagator described in section 4 
of paper I where the ghost particle self-energy is written as Vl{k) = k'^Cl{k'^) and combined with which comes from 
the kinetic energy operator □. Therefore, in the general gauge, the ghost particle kinetic energy term in Eq. (3.10) 
may be extended to the following form by directly adding the ghost particle mass term 

Y^^YcC^nC" + CV^C" (3.22) 

where the factors Y^^ and Yc arise respectively from the change of the variables (7° and C" which are caused by the 
counterterm contained in the second term in Eq. (2.27), as described in Eqs. (3.3) and (3.4). With the definitions 
given in Eqs. (3.12) and (3.15), Eq. (3.22) can be represented as 

Z3(7»(n + Zg- V')(^" = (° + Mo)Co (3-23) 

where 



/i = VZ3Mo (3.24) 

which just is the relation given by Eq. (4.30) in paper I. Clearly, in the general gauge, the ghost particle kinetic energy 
terms in Eqs. (3.13) and (3.17) should be replaced by the terms on the LHS and RHS of Eq. (3.23), respectively. 

IV. GLUON THREE-LINE VERTEX 

In the preceding section, the S-T identity shown in Eq. (3.14) for the renormalization constants was derived 
from the general proof of the renormalizability of the QCD with massive gluons. Originally, the same identity given 
in the QCD with massless gluons was derived from the W-T identities satisfied by the proper vertices and their 
renormalizationt^^'^^1. Following the similar line, we will re-derive the identity in Eq. (3.14) so as to give a check of 
the identity and its derivation. Since the difference between the massive QCD and the massless QCD mainly lies in 
the gluon field, in this section and the next section, we limit ourself to take the gluon three-line and four-line proper 
vertices as examples to show how to derive W-T identities satisfied by the vertices and discuss their renormalization. 
To derive W-T identities for the vertices, it is suitable to employ the identities written in Eqs. (2.14) and (2.15) 
though, we would like to start from the identities satisfied by the generating functional of Green functions in order 
to see, by the way, the difference between the Green functions for the massive gluons and the massless ones. Let us 
begin with derivation of an identity satisfied by the gluon three-point Green function from the following W-T identity 

1^^ 6Z[J] , [ ^ 5-'Z[J,K,u] , 

-^x ,-r„„,r^, + j dyJ" {y) ,r^„,r^,,^,^„r.^ K=u=0 = (4.1) 



a ^^J«''(a;) J ^ ' 5K'^{x)6vP''{yy 



and the ghost equation 



■p,x S^Z[J.K.u] I ■ 2 5^Z[J.K,K] | _ 

^"fi Su^(x)SK''(y) \K=u=0 + «M SK'^{x)SK''{y)\K=K=0 j-^^) 



+S'''>5^{x -y)Z[J]=0 
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which were derived in section 4 of paper I ( Hereafter, /i denotes the unrenormaHzed mass of ghost particle and fiji 
designates the renormahzed mass of the particle). By taking successive differentiations of Eq. (4.1) with respect to 
the sources J^{y) and J^(z) and then setting the sources to vanish, one may obtain the W-T identity obeyed by the 
gluon three-point Green function which is written in the operator form as 

a 

+ < Q+\T*[C-{x)Al{y)Dl^{z)C''{z)]\0- > (4.3) 

where 

G%ix,y,z) =< 0+\T[A;{x)At{y)Al{z)]\0- > (4.4) 

is the three-point Green function mentioned above. The identity in Eq. (4.3) will be simplified by a ghost equation 
which may be derived by differentiating Eq. (4.2) with respect to the source Jxiz) 

dH < 0+\T*{D';f{x)C''{x)C'^iy)Aliz)}\0- > 
< 0+\T[C''{x)C'>{y)Al{z)]\0-] >= 0. 

Taking derivatives of Eq. (4.3) with respect to y and z and employing Eq. (4.5), we get 

dlid;d^G;i\{x, y, z) = aii^dl^'^liz, X, y) + a,^G^^(y, x, z)} (4.6) 

where 

G'''';{x,y,z) =< 0+|T{C"^(x)(7''(y)i^(2)}|0- > . (4.7) 

In the Landau gauge or the zero-mass limit (/x = 0), Eq. (4.6) reduces to 

«5.^G«t'x(a=,2/,^) = (4.8) 

which shows the transversity of the Green function. From Eq. (4.6), we may derive a W-T identity for the gluon 
three-line vertex. For this purpose, it is necessary to use the following one-particle-irreducible decompositions of the 

Green functionsl^^'""^^! 

Gt\{x, y,z) = J d^x'd^y'd^z'iDl^,{x - x') 

xiDlUy-yyDTyi^-^'K'':'c''i^',y',z',) ^ ' ^ 

and 

G°'''l(x,y,z) = J d'^x'd'^y'd'^z'iA""' (x - x')^''''''' {x' ,y' , z') , . 

xiM>'\y' - y)iDi,l{z' - z) ^ ' ^' 

where il?"^", (.t — .t') and iA"" [x — x') are the propagators for the gluon and the ghost particle respectively which were 

derived in paper I, ^abci^-' V' ^) ^'^'^ r"''^(x, y, z) are the three-line gluon proper vertex and the three-line ghost-gluon 
proper vertex respectively. They are defined as'^^'^^l 

and 

"""^-'^'-) = ^.C-(x).CH.)Mc^(.) l--° ^'-''^ 

hereJ stands for all the external sources. Substituting Eqs. (4.9) and (4.10) into Eq. (4.6) and transforming Eq. (4.6) 
into the momentum space, one can derive an identity which establishes the relation between the longitudinal part of 
three-line gluon vertex and the three-line ghost-gluon vertex as follows 
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p^q^k^Af-^ip, q, k) = -!^x{P^)[x(k'')q''A--t{k,p, q) ,^ ^3. 

+X(g2)fc^A''-^-(g,p,A:)] 



where we have defined 



and 



x(p') = + nL(p')] - + + - + (4.15) 

here IIl(p^) and ^^(p^) are the self-energies appearing in the longitudinal part of gluon propagator and the ghost 
particle propagator (see Eqs. (4.11) and (4.22) in paper I). 
Obviously, in the Landau gauge, Eq. (4.13) reduces to 

p''(Z''fcXtA(P,9,fc) = (4.16) 

which implies that the vertex is transverse in this case. In the lowest order approximation, owing to 

X(p') = 1 (4.17) 

and 

A(°)'"';(p,g,A;)=5r''=p^, (4.18) 
the RHS of Eq. (4.13) vanishes, therefore, we have 

p''g-A;^A(°);:tA(P>9,fc) = 0. (4.19) 

This result is consistent with that for the bare three-line gluon vertex given by the Feynman rule. 

Now, let us discuss renormalization of the three-line gluon vertex. At first, we note that the notation used here 
is different from the previous section. In this section, a renornalized quantity will be marked with a subscript R\ a 
quantity without the subscript R represents the unrcnormalizcd one. From Eq. (3.15) or the renormalization of the 
gluon and ghost particle propagators as described in paper I, wc see 



A-^^ix) = (x). 

According to these relations and the definitions given in Eqs. (4.11), (4.12) and (4.14), we find 

K^\ip. Q, k) = Z-^/^A^;ll{p, q, k), 
A-^lip, q, k) = Z^^Z-^'^K^-\{p, q, k). 



(4.20) 



(4.21) 



Applying these relations, the renormalized version of the identity written in Eq. (4.13) will be 

P^'fk'^^^fAiP^'l^k) = -^XR{p')[XR{k^)q''A^'t{k,p,q) .4 22) 

+Xfl(9')fc^A/«jf(g,p,fc)] 

where an = Z^^a, = \fz-j,\x and 

is the renormalized expression of the function x(fc^) which may be obtained by substituting the relation (see Eq. 
(4.23) in paper I) 

n^(fc2) = i^^^^k^) (424) 
^ ' fc2[l + 0(fc2)] ^ ' 
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into Eq. (4.15) and then using the definitions 

= + n%k'^), = 1 + Qnik'^) = ZsQ^k'^) (4.25) 

which were ever shown in section 4 of paper I. At the renormahzation point chosen to be — = k"^ = fij^, we see, 
XBil^ji) = 1- In this case, the rcnormalized ghost-gluon vertex takes the form of the bare vertex so that the RHS of 
Eq. (4.22) vanishes, therefore, we have 

p^q'^k^Aj.^^tlip, q, A;)|p.=,.=fe.=^2^ = 0. (4.26) 

Ordinarily, one is interested in discussing the renormahzation of such three-Une vertices that they arc defined from 
the vertices defined in Eqs. (4.11) and (4.12) by extracting a couphng constant g. These vertices are denoted by 
A^^'3^(p, q, k) and h°'^^{p, q, k). Commonly, they are rcnormalized in such a fashionl^-^'-^^l. 

AfA(p, q, k) = Zj'I^b^^.ip, q, k), 
K-\{p, q, k) = Z^^K-j^^ip, q,k) 



where Zi and Zi are referred to as the renormahzation constants for the three-line gluon vertex and the ghost-gluon 

T ab< 



vertex, respectively. It is clear that the W-T identity shown in Eq. (4.13) also holds for the vertices A"^'3^(p, g, fc) 



and A"''^(p, q, k). So, when the vertices A^^'3^(p, q, k) and h.'^\{p, q, k) in Eqs. (4.13) are replaced by AJ^^aCP' 9' ^) ^"^^ 
K°-\{p, q, k) respectively and then Eq. (4.27) is inserted to such an identity, we obtain a renormalized version of the 
identity as follows 

P>^q''k^l4l\{p,q,k) = -|i|^x«(p2)[^^(fc2) 
xq'^K^<^^{k,p, q) + XK(9')fc^A^f<\-(9,p, k)]. 

When multiplying the both sides of Eq. (4.28) with a renormalized coupling constant qr and absorbing it in the 
vertices, noticing 



^t\{p,q,k)=gRKf\{p,q,k), 



(4.29) 
we have 



p^q-'k^A^;ll{p,q,k) = -||^Xfl(p')[XK(fc') .4 30) 

X9'^A^-^(fc,p, (?) + xk((Z^)A;^A/<\^(9,P, fc)]. 

In comparison of Eq. (4.30) with Eq. (4.22), we see, except for the factor ZiZzZ^^Z^^, the both identities are 
identical to each other. Prom this observation, we deduce 

I = I 

This is the S-T identity which coincides with the one given in Eq. (3.14). 

V. GLUON FOUR-LINE VERTEX 

By the similar procedure as deriving Eqs. (4.3) and (4.4), the W-T identity obeyed by the gluon four-point Green 
function may be derived by differentiating Eq. (4.1) with respect to the sources J^{y), J\{z) and J^{u). The result is 
as follows^ 

ia^fG«tf.(x,y,z,«) 

=< 0+\T*[C-{x)Dl^iy)C%y)Al{z)Ai{u)]\0- > 

+ < 0+\T*[C^{x)Al{y)Dl%z)C%z)Ai{u)]\0- > ^ ' ' 

+ < 0+\T*[C''{x)A'l{y)Al{z)D^^{u)C%u)]\0- > 
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where 

Gf,^t{x,y,z,u) =< 0+\T[A^^{x)Ai{y)Al{z)Ai{u)]\0- > (5.2) 

is the gluon four-point Green function. The accompanying ghost equation may be obtained by differentiating Eq. 
(4.2) with respect to the sources Jxi^) and J!^{u). The result is 

< 0+\T*[D';f{x)C''ix)C\y)A 
^^2Qabcd^^^ z, u) = -S'^H^ix - y)Dl''^{z - u) 



< 0+\T*[D-^ix)C%x)Cbiy)Aliz)Aiiu)]\0- > 3^ 



where 

G''\i{x,y,z,u) =< 0+\T[C'^{x)d\y)Aliz)Aiiu)]\0- > (5.4) 

is the four-point gluon-ghost particle Green function. Differentiation of Eq. (5.1) with respect to the coordinates y, 
z and u and use of Eq. (5.3) lead to 

hd^dp^dlGH^^^ix, y, z, u) = S'^^S^x ~ y)d^d:Dli{z - u) 
+<5«^(54(x - z)d^dlDti{y -u)+ S'^'^S^x - u)d-yd^Dl\{y - z) 

+fi^{d^dlG>"'^xiiy, x, z, u) + d-ydlG'^'^tiiz, x,y,u) ^ ' ' 

+d^d^G'^''^l{u,x,y,z)}. 

It is noted that the four-point Green functions appearing in the above equations arc unconnected. Their decompo- 
sitions to connected Green functions are not difficult to be found by making use of the relation between the generating 
functionals Z for the full Green functions and W for the connected Green functions, Z = exp{iW). The result is 

G'Jlfri^' ^' ^) = Gliti^, y, z, u)e - D-U^ - y)Dli{z - u) 

-Dllix - z)Dli{y -u)- D-iix - u)Dll{y - z) ^^•''^ 

and 

G'^'liix, y, z, u) = G'^'^xiix, y, z, n)e - A-\x - y)D'ii{z - u). (5.7) 

The first terms marked by the subscript "c" in Eqs. ( 5.6) and (5.7) arc connected Green functions. When inserting 
Eqs. (5.6) and (5.7) into Eq. (5.5) and using the W-T identity satisfied by the gluon propagator which was derived 
in section 4 of paper I 

dUd^Df^ix -y)- aa^A'^^x - y) = -a8''''5\x - y), (5.8) 



one may find 



d^d-dzdlGf^'t y. z, u\ = ati^dp^G^-^^xiu, x, y, z), 

+d^dlG-'^^f{z,x,y,u)e + d^d:,G'>--xi{y,x,z,u)c}. ^ ' ' 

This is the W-T identity satisfied by the connected four-point Green functions. In the Landau gauge, we have 

dlid^d^dlGUti^, y, z, u)e = (5.10) 

which shows the transversity of the Green function. 

The W-T identity for the four-line gluon proper vertex may be derived from Eq. (5.9) with the help of the 
following one-particle -irreducible decompositions of the connected Green functions which may be found by the 
standard procedure'^^'^^l, 

Gl'l''^^ixi,X2,X3,Xi)c 



= 111 d'yiDZ,{x, - yi)Dit ix2 - 2/2)r^:r'X''(2/i.y2,y3,y4) 



xDl,l{ys-Xs)D^,,^,iyi-x^) 

+iJU dSid'zM'i.ixr - yMU^2 - y2K'/,::{yi,y2,ys) 



fj.fl - 

' , , (5-11) 

x7?- (yg - ^i)r^,e.; {z,,Z2, z:i)Dl,liz2 - x^)Di,i{z, - X4) 

(a;i - y^)Dli (x^ - y2)r'::^::{yi , ^2, ?/3)i?^^: (ya - ^i) 

Xr^:^;j;(^l, ^2, ^3)A''t(^2 - X2)Df,iiz3 - X^) 

+dZ{x2 - yi)Dli{x3 - y2)T'(,^,:{yuy2,yz)Dli,{y^ - z^) 



13 



and 

= -i/ n rf'yiA""'(a;i - yi)T^,,,^;s:{yuy2,y3,y4)^'''Hy2 - X2) 

xDi,l{ys-xs)Df4{y4-X4) 

+iJUdVzi{A--'{xi - yi)r„,/;,(yi,y2,y3)A-'(j/2 - ^i) 

xDf4{ys - X4)T^,,,^;{z,,Z2, z^)A^'\z2 - X2)Di,\{z^ - x^) 
+A-'(a;i - ?/i)r^,i;(j/i,2/2,y3)A-'(y2 - zi)Di,\{y:, - xs) 

Xr,,(,7-;,(Z1,Z2,Z3)A'''''(02 - X2)Df4{z3 - X4) 

+A--'{xi - yi)T^,^^^{yi,y2,y^)A!''\y2 - X2)D;<;;,{ys - ^i) 

>^K'c'l'izi,Z2, Zs)Di,l{z2 - Xs)Df'4{zs - X4)} 

where r'jfj'^^{xi,X2,X3,X4) is the four-Une gluon proper vertex and T"'''^f{xi,X2,X3,X4) is the four-line ghost-gluon 
proper vertex. They are defined as 



r 



^lt'AT(^l'^2,a;3,a;4) — i sA^i'(xi)5Al>"(x2)SA''^(x3)SA''^(xi)\j=0^ (K-io\ 
1 AtV-^1''^2,^3, ^-4; — iSC''{xi)5C"'{x2)SA''>-{x3)5Ad^{x4)^-^=^' 



When substituting Eqs. (5.11) and (5.12) into Eq. (5.9) and transforming Eq. (5.9) into the momentum space, one 
can find the following identity satisfied by the gluon four-line proper vertex 



c d 

'4 



k^k^,k^klA;itiki,k2,h,h) = ^^l^^ I 

a c d h \ ^ I a d b c 
ki k3 ki k2 ) \ ki /c4 k2 k^ 



(5.14) 



where 



abed 
ki k2 ks k4 

■ik'(k}iK%{kr,k2, -ih + A;2))^?:;(fci + k2)k^klAlli{-{k3 + ki), ki, ki) 

+ '^x{kl)x{kl)[ik^klk^-li{k2M,hM) (5.15) 

-A^«^(fe, fci, -(fci + fe))i?:;(fcl + fe)fc3^fcJAj^^(-(fc3 + ki), k3,ki) 

~klh'-i{k2,-{k2 + ki), ki)A-f{k2 + ki)k^A^'^^i-iki + ks), fci, fc3) 

-k^A''^l{k2, -{k2 + ks), fc3)A^/(fc2 + fc3)fcIA/«^(-(fci + ki), fci, ki)]. 

The second and third terms in Eq .(5.14) can be written out from Eq. (5.15) through cyclic permutations. In the 
above, we have defined 

r«tt (fci, fc2, ks, ki) = i27r)^S\j:U h)Kl'^riki, k2, ks, ki), 

T-\i{kuk2, ks, ki) = {lTT)H\Y.t=i ki)A-\i{k^, k2, ks, ki). ^ ' ' 

In the lowest order approximation, we have checked that except for the first term in Eq. (5.15) which was encountered 
in the massless theory, the remaining mass-dependent terms are cancelled out with the corresponding terms contained 
in the second and third terms in Eq. (5.14). Therefore, the identity in Eq. (5.14) leads to a result in the lowest order 
approximation which is consistent with the Feynman rule. 

The renormalization of the four-line vertices is similar to that for the three-line vertices. From the definitions given 
in Eqs. (5.13), (5.16) and (4.20), it is clearly seen that the four-line vertices should be renormalized in such a manner 

Ktiki,k2, ks, ki) = Z^^A^;txr{ki,k2, ks, ki), 
A-Yriki,k2, ks, ki) = Zs'Z^^Aj,-t'l^^ikuk2, ks, ki). ^ ' ' 

On inserting these relations into Eqs. (5.14) and (5.15), one can obtain a renormalized identity similar to Eq. (4.22), 
that is 
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a c d ^ \ _|_ f d b (■ 
ki ks ki k2 I ^ V ^1 ^4 A;3 



(5.18) 



where 



a 



bed 
ki k2 k^ ki 

-<k:^^^t^laikl,k2, -(fcl + k2))D^Pj.{k^ + k2)k^klAif^^{-{k3 + k4),k3,k4) 

+ '-^XR{kl)xR{kl)[ik^klk'^^-\^{k2MMM) (5.19) 
-A'^'^,{k2M, -{ki + fc2))£'fle/(fci + k2)k^klAl;l{-{ks + k^), ks, k^) 

-klkYAk2, ~(fc2 + fc4), fc4)A^^(fc2 + ki)k^kl,"\{-{ki + ks), fci, ks) 
-k^A'g\{k2, -{k2 + ks), ks)^ri {k2 + k3)klA^^^ii-{ki + ki), ki,ki)] 

Wc can also define the vertices A'^J'^^{ki, k2, k^, k^) and A°'^'^^^{ki,k2,ks,k4) by taking out the coupling constant 
squared from the vertices A.'^^^^{k\,k2,k3,ki) and A°-^^{k\,k2,k3,ki), respectively. The renormalization of these 
vertices are usually defined by^'^ 

A;tA.(fci, k2, ks, ki) = Z^'Aj,-Xik^, k2, ks, ki) 

A-YAki,k2, ks, ki) - Z^'Aj,-'^Ak,,k2,ks, ki) ^ ' ^ 

where Zi and Zi are the renormalization constants of the four-line gluon and ghost-gluon vertices respectively. 
Obviously, the identity in Eqs. (5.14) and (5.15) remains formally unchanged if we replace all the vertices Aj in the 
identity with the ones . Substituting Eqs. (4.27) and (5.20) and the following renormalization relations for the 
gluon and ghost particle propagators (which were given in section 4 of paper I) 

iA''\k) = ZaiAf (fc) 



iD;i{k) = ZsiD^;t{k) (5.21) 
into such an identity, one may write a renormalized identity similar to Eq. (4.28), that is 



where 



k'^k!^klklAg^;iUki,k2,ks,ki) = ^ni ^ ^ 



^ , a c d b \ I a d b c 



ki ks ki k2 J \ ki ki k2 ks 



abed 

ki k2 ks ki 

= ^{-ik'tkfA-^;^^{k,,k2, -(fci + k2))D^:..{k, + k2)klklAl:l{-{ks + ki), ks, ki)} 



(5.22) 



+ '-ixR{kl)XR{kl){i^k^kfA}t^t{k2,ki,ks,ki) ^523) 

-|rl:^i?"'-(^2, -(fci + k2))D^p^^.{k, + k2)k^kiAl;l^{~{ks + ^4), ^3, ki) 

-^M~^Tr{k2,~{k2 + ki),ki)/A^i{k2 + ki)kfAl''\{-{k^ + ks),ki,ks) 
+k^V>^\{k2,-{k2 + ks),ks)Arl{k2 + ks)klV^''^{-{ki + ki),ki,ki)\} 

Multiplying the both sides of Eqs. (5.22) and (5.23) by g\, according to the relations given in Eqs. (4.29) and in the 
following 

^ZtxrikiM. kz, ki) = glAj,f,tiki, k2,ks,ki) 

^R'irikl,k2, ks, ki) = glAj,-\^{ki,k2,ks, ki) 
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we have an identity which is of the same form as the identity in Eqs. (5.22) and (5.23) except that the vertices in 
Eqs. (5.22) and (5.23) are all replaced by the vertices A^. Comparing this identity with that written in Eqs. (5.18) 
and (5.19), one may find 

^ = 1,^ = 1,^ = 1,^ = 1 (5.25) 
Zl ' Z^Zi ' Z^Zr ' z^zl ^ ' 

which lead to 

^ = £l = ^,£i = ^ = S. (5.26) 

Zs Zo, Zy Z4 

This just is the S-T identity which is consistent with that given in Eq. (3.14). 



VI. EFFECTIVE COUPLING CONSTANT 



To concretely demonstrate the renormalizability of the massive gaiige field theory, in this section, we take one- 
loop renormalization of the QCD with massive gluons as an example. The renormalization is carried out by the 
renormalization group approach^^^"^^! and gives an exact one-loop effective coupling constant of the QCD. As argued in 
our previous paperl^^"^^!, when the renormalization is carried out in the mass-dependent momentum space subtraction 
schemeP*~^^l, the solutions to the renormalization group equations (RGE) satisfied by renormalized wave functions, 
propagators and vertices can be uniquely determined by boundary conditions of the renormalized wave functions, 
propagators and vertices. In this case, an exact S-matrix element can be written in the form as given in the tree- 
diagram approximation provided that the coupling constant and particle masses in the matrix element are replaced 
by their effective (running) ones which are given by solving their renormalization group equations. Therefore, the 
task of renormalization is reduced to find the solutions of the RGEs for the renormalized coupling constant and 
particle masses. Suppose Fji is a renormalized quantity. In the multiplicative renormalization, it is related to the 
unrenormalized one F in such a way 

F = ZfFr (6.1) 

where Zp is the renormalization constant of F. The Zp and Fji are all functions of the renormalization point /i = /Uqc* 
where is a fixed renormalization point corresponding the zero value of the group parameter t. Differentiating Eq. 
(6.1) with respect to \x and noticing that the F is independent of /x, we immediately obtain a renormalization group 
equation (RGE) satisfied by the function F^~'^'^ 

pi^+-iFFR = Q (6.2) 
ail 

where is the anomalous dimension defined by 

7F = At^ln.^F- (6.3) 

Since the renormalization constant is dimensionless, the anomalous dimension can only depend on the ratio /3 = — 

where m/j denotes a, renormalized mass and 7^= TfI.?/?'/^) which g/j is the renormalized coupling constant and 
depends on \i. Since the renormalization point is a momentum taken to subtract the divergence, we may set \i = ji^X 
where A = e* which will be taken to be the same as in the scaling transformation of momentum p = poX. In the 
above, /io and po are the fixed renormalization point and momentum respectively. When we set F to be the coupling 
constant g and noticing = A^, one can write from Eq. (6.2) the RGE for the renormalized coupling constant 

A^^+7.(%k(A)=0 (6.4) 



with 



73=M^ln^g. (6.5) 
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The renormaHzation constant Zg is commonly defined by (see Eq. (3.15)) 



= = ^ (6-6) 
Z^Zi 



where the last equality is given by using the identity in Eq. (4.31). Here we would like to choose the expression of Zg 
given by the last equality in Eq. (6.6) to evaluate the anomalous dimension 7g(A). As denoted in Eqs. (4.25) in paper 
I and Eq. (4.27), the renormalization constants Z3, Z3 and Zi are determined by the gluon self-energy, the ghost 
article self-energy and the ghost vertex correction, respectively. At one-loop level, the gluon self-energy is depicted in 
Figs, (la)-(ld), the ghost article sclf-cncrgy is shown in Fig. (2) and the ghost vertex correction is represented in Figs. 
(3a) and (3b). According to the Feynman rules which are the same as those for the massless QCD1^°1 except that 
the gluon propagator and the ghost particle one are now given in Eqs. (4.14) and (4.10) in paper I, the expressions 
of the self-energies and the vertex correction are easily written out. For the gluon one-loop self-energy denoted by 
— zn^^(A;), one can write 

4 

where n|^^j°^(A:) — n|^'|]°''(fc) represent the self-energies given in turn by Figs.(la)-(ld). They are separately represented 
in the following: 

n|ij (fc) = Z(5"''|g2 / [j2_„2+j£][(;+ft)2_„2+j£] [gp.\{l + 2fc)p - g\p{2l + k)^ /g gx 

iV [k)- Id 35 y ^2;r)4[(« + /fc)2-m2 + i£][F-m2 + i£]' ^^"^^ 

nifj"^(fc) = -iS^%' J (2^^^ Q2_C + i,fr '^9^P - 9,P9X.) (6.10) 

and 

T-r(4)a6/7\ -sabl 2 f d^l 1 

yyp.v [Kj--10 25 J [(l-ky^-M'^+is][l'^-M'^+ie] (Qll) 

xTr[7^(l-k + M)7,(l + M)] ^ ' ^ 

where 1 =7^;a, k =7^fcA. In the above, f'^'^^jbcd ^ g^ab ^^^^ Tr{T''T^) = ^d"^ have been considered. It should 
be noted that in writing Eqs. (6.8)-(6.10), wc choose to work in the Feynman gauge for simplicity. This choice is 
based on the fact that the massive QCD, as proved in paper I, is an unitary theory, that is to say, the S-matrix 
elements evaluated from the massive QCD are independent of gauge parameter. Therefore, we are allowed to choose 
a convenient gauge in the calculation. From Eqs. (6.8)-(6.11), it is clearly seen that 

4 

n;t(fc) = ^''X^(fc) = (6.12) 

i=l 

By the dimensional regularization approach^^J, the divergent integrals over / in Eqs. (6.8)-(6.11) can be regularized 
in a n-dimensional space and easily calculated. The results are 

n^jj(A:) = -| (^j^ /o s[kix{x-i)+m-^Y {9p.i'[^'^x{x - 1) ,g 
+b)k'^ + 9rn^]+2[bx{x-l)-l\k^K}, 

(fc) = I (^;^ /q e[k'^x(x-l)+m'^Y {[^^^(^ ~ -*") (6.14) 

+"^^]5/il/ + 2a;(a; - l)/;;^fc^}. 
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and 

where £ = 2 — ^ ^ when n 4. In Eqs. (6.13)-(6.16), except for the e in the factor l/e[k'^x{x — 1) + m'^Y 
l/e[k'^x{x — 1) + M'^Y, we have set £ — > in the other factors and terms by the consideration that this operation does 
not affect the calculated result of the anomalous dimension. According to the decomposition shown in Eqs. (4.15) 
and (4.16) in paper I and noticing = P^" + P^" where Pl^" = {g>"' - ^) and P^" = it is easy to get the 

transverse part of njui/(fc) from Eqs. (6.13)-(6.16) and furthermore, based on the decomposition denoted in Eq. (4.20) 
in paper I, i.e., nT(fc^) = A;^ni(A;^) +m^n2(A;^), the functions ni(fc^) and n2(A;^) can be written out. The results are 

n, ik') - t dx\ 15[2x(x--l) + l] ^ ^x{x - 1) 

{^-Kf }^ ^2e\kH(x-\) + m^Y^ y^e{k''x{x-\) + MlY^ 



and 



g2 /-i 27 

n2(fc2) = -yf-TTjl / dx-rr^, 7T ^ - - — 7}. (6.18) 

^ ' {^ttY Jo e[k^x{x-l) + m^Y 



It is clear that the both functions ni(fc^) and Il2{k^) are divergent in the four-dimensional space-time. When the 
divergences are subtracted in the mass-dependent momentum space subtraction schemet^^"^"'^!, in accordance with 
the definition in Eq. (4.25) in paper I, we immediately obtain from the expression in Eq. (6.17) the one-loop 
renormalization constant Z3 as follows 

^3 = 1-01(^2) 

,^ (6.19) 



g\ fl , f 15[2x(x-l) + l] ixjx-l) 1 

i=l * 



Next, we turn to the ghost particle one-loop self-energy denoted by —iQ,"'^{q). Prom Fig. (2), in Feynman gauge, 
one can write 

0«^(g) = I \^'^7,ll (6.20) 

^ J (2^)^ [{q - l)^ - + ie][P - + le] ^ ' 

By the dimensional regularization, it is easy to get 
where 



l_ 3(.-l) 

(47r)2yo e[q^x{x-l)+m''Y 



^(^') - I dx ^^^jr o^. - (6.22) 



According to the definition given in Eq. (4.25) in paper I and the above expression , the one-loop renormalization 
constant of ghost particle propagator is of the form 

= 1 - Cli,^) = 1 - (1)2 '^^ ei.^xit'lHm^Y - ^'-''^ 

Now, let us discuss the ghost vertex renormalization. In the one-loop approximation, the vertex defined by 
extracting out a coupling constant is expressed as 

Af ^(P, q) - r"'px + Ati^ip, q) + Kt^ip, q) (6.24) 

where the first term is the bare vertex, the second and the third terms stand for the one-loop vertex corrections shown 
in Figs. (3a) and (3b) respectively. In the Feynman gauge, the vertex corrections are expressed as 
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ia6c/^„x_ ,fabc^ 2 f P ' {q - 1){P - I) \ 



and 



lAk-f^.y; J 2^ J {2Tryp-m^ + ie][{p-l)^-m^+ie][{q-l)^-m'^ + ie] ^ ' 



Aabc/ X ^ ,-fa6c3 2 /" ? I ■ (p - Q - l)p\ - p ■ Iqx + p ■ i2q - p + l)lx 

^^2X\P,q) 'J 2^ J {2tt)^P - m'' + ie][ip - l)^ - + ie][iq - l)^ - m-' + ie] ^ ' 

where f'^'^d.jebfjdfc _ _ 3j:abc j^g^g ^gg^ notcd. By employing the dimensional regularization to compute the above 
integrals, it is not difficult to get 

Afffe») = - ^/^^'M + ..B.M] - i«) (6.27) 



where 



and 



where 



^xy = P^xy{xy - 1) + q^[{x - ify + (x - l)]y - 2p • qx{x - l)y'^ + m^, 

Ai{p, q)^{p- q[l + {x- l)y] - p^xy}{l - xy)y, (6.28) 
Bi {p, q) = {p- g[l + (x - l)y] - p^xy}{x - l)y^ 



Kx{P,q) = r''\7fyi fdx ['dy{f^ + -^\pxA2{p,q)+qxB2{p,q)]-lpx} (6.29) 



-42(p,?) = {p'^{2xy - x2y2 „ 1) _ q-2[(^x _ 1)^ _ l](a; - l)y 

+p-q[2-{3x-2)y + 2x{x-l)y^]}y, (6.30) 
-B2(p, q) = \p-q{x-l)- p'^x]y'^. 

The divergences in the both vertices Al\'^{p,q) and A2^'^(p, g) may be subtracted at the renormalization point p"^ = 
q^ = fj? which implies k ~ p — q = being consistent with the momentum conservation held at the vertices. Upon 
substituting Eqs. (6.27) and (6.29) in Eq. (6.24), at the renormahzation point, one can get 

If^ip, q) 1^2=,.=^.= r^'^pxil + Zi) = V^'^A (6.31) 

where 

{^-Kf Jo ^^e\pi^x{x-l) + 'm?Y ^J:'x{x-l)+m? 4^ ' 

which is the one-loop renormalization constant of the ghost vertex. 

Now we are ready to calculate the anomalous dimension 7g(A). Substituting the expressions in Eqs. (6.6), (6.19), 
(6.23) and (6.32) into Eq. (6.5), it is easy to find an analytical expression of the anomalous dimension 7g(A). When 
we set ^ = ^ and ^ = ^ with defining P = and Pi = ^ (here we have set /xq = A), the expression of 7g(A), in 
the approximation of order g^, is given by 

7,(A) = lim[.AlnZ, - ,|-lnZ3 - l.^lnZs] = ^F(A) (6.33) 

where 

A+v^A^^ 3A.^l^+A= X^^X?^, X+^>^^^ 

in which iVy denotes the number of quark flavors. We would like to note that the flxcd renormalization point A in /? 
and Pi can be taken arbitrarily. For example, the A may be chosen to be the mass of the quark of N f -th flavor. In 
this case, /? = rajM^j and pi = Mi/M^^. In practice, the A will be treated as a scaling parameter of renormalization. 
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With the 7g(A) given above, the equation in Eq. (6.4) can be solved to give the effective coupUng constant as 

follows 

= T^g^ (6.35) 

where an = Q!k(1) and 

G(A) = / - = M>^) - M^) - 3 E[^2(A) - ^^(1)] (6.36) 

in which 



^i(A) = [(19 - + VA^-4/3^ J A-VA^-4/3^ + ^' ^'-''^ 



.UA) = (f.f)V^ln^±^-^ 

and <^i(l) = <pi(A) |a=i, ¥'2(1) = '^2(A) Ia=i- In the large momentum limit (A 00), we have 

G(A) = (11- ^Ar/)lnA. (6.39) 

This just is the result for massless QCD which was obtained previously in the minimal subtraction scheme[^*~^°l. It 
should be noted that the expressions in Eqs. (6.34), (6.37) and (6.38) are obtained at the timelike subtraction point 
where the A is a real variable. Wc may also take spacclikc momentum subtraction. For this kind of subtraction, 
corresponding to /U ^ i/j., the variable A in Eqs. (6.34), (6.37) and (6.38) should be replaced by iX where A is still a 
real variable. It is easy to see that the function in Eq. (6.39) is the same for the both subtractions. 

The behavior of the function a/j(A) is graphically described in Figs. (4) and (5). Figs. (4) and (5) represent 
respectively the effective coupling constants obtained at the timelike subtraction point and the spacelike subtraction 
point, where we take the flavor Nj = 3 as an illustration. The parameters aji, rn and A are taken to be a^. = 0.2, 
m = 600Mey, A = SOOMeV. The masses of up, down and strange quarks are taken to be constituent quark ones, 
i.e., M„ = Met = 350Mey and Mg = bOOMeV. For comparison, we show in each of the figures three effective coupling 
constants which are obtained by the the massive QCD, the massless QCD and the minimal subtraction, respectively. 
These effective coupling constants are respectively represented by the solid, dashed and dotted lines in the figures. 
Prom Fig. (4) it is seen that the effective coupling constant given by the massive QCD is an analytical function with 
a maximum at A = 1.346, the effective coupling constant given by the massless QCD is also an analytical function 
with a peak around A = 1, but the effective coupling constant given by the minimal subtraction has a singularity 
at A = 0.1746. Fig. (5) indicates that the effective coupling constant given by the massive QCD, similar to the one 
given by the minimal subtraction, has a Landau pole at A = 0.1845 which implies that the coupling constant is not 
applicable in the region A < 0.1845. However, if the gluon mass is taken to be m < 425.75Mey, we find, the Landau 
pole disappears and the effective coupling constant, analogous to the effective coupling constant given by the massless 
QCD , becomes a smooth function in the whole region of momentum as illustrated by the dotted-dashed line in Fig. 
(5) which represents the effective coupling constant given by taking m = 425.75MeV^. 

In this paper, we limit ourself to show the derivation and result of the effective coupling constant only. The 
effective gluon mass and quark mass can be derived in the similiar way. All these effective quantities, as expected, 
are of asymptotic behaviors. 



VII. CONCLUDING REMARKS 



The derivation and the result stated in section 2 clearly show that the divergences appearing in the perturbative 
calculations for the massive non-Abelian gauge field theory can indeed be eliminated by introducing a finite number of 
counterterms as shown in Eqs. (2.27)-(2.31) and (2.34). Saying cquivalcntly, as shown in section 3, these divergences 
may be absorbed into a finite number of renormalization constants and wave functions and thus be removed by 
redefining the wave functions and the physical parameters. In view of this, according to the general argument of 
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renormalizability, wc can say that the renormahzabihty of the QCD with massive gluons is absolutely no problem. 
The renormalizability of the massive QCD was illustrated in section 6 where we have derived a rigorous one-loop 
effective coupling constant and shown that the corresponding result obtained in the previous literature[^*~^°l is 
only an approximate one given in the large momentum limit. Thus, contrary to the prevailing concept that it is 
impossible to build a renormalizable massive gauge field theory without recourse to the Higgs mechanism I^-io), we 
have succeeded in establishing such a theory which is renormalizable. The basic idea to achieve this success is the 
consideration that the massive gauge field must be viewed as a constrained system in the whole space of the full vector 
potential. Therefore, to construct a correct quantum theory for such a system, the unphysical degrees of freedom 
contained in the massive Yang-Mills Lagrangian must be eliminated by introducing appropriate constraint conditions 
on the gauge field and the gauge group. In doing this, the requirement of gauge-invariance must be respected in the 
whole process of the construction of the theory. Particularly, in the physical space defined by the Lorentz condition, 
only infinitesimal gauge transformations are necessary to be considered. These essential points were not realized 
clearly and handled correctly in the previous studies. In the earlier works of investigating the massive non-Abclian 
gauge field thcoryl^^^^l, as mentioned in Introduction, authors all started with the massive Yang-Mills Lagrangian 
and considered that this Lagrangian itself forms a complete description of the massive gauge field dynamics. When 
using this Lagrangian to construct the quantum theory, they found that except for the neutral vector meson field 
in interaction with a conserved current, the theory is nonrenormalizable because of the presence of the mass term. 
However, as pointed out in Ref.[l], when the Lorentz condition is introduced and thereby the infinitesimal gauge 
transformations are taken into account only, the problem of nonrenormalizability will all disappear in those works. 

Another point wc would like to emphasize is that for the massive gauge field theory, the ghost particle acquires 
a spurious mass /U in general gauges. This mass term is necessary to be introduced so as to compensate the gauge- 
non-invariance of the gauge boson mass term and preserve the effective action to be gauge-invariant. As we see, the 
occurrence of this mass term in the theory is essential to guarantee the theory to be self-consistent. Otherwise, for 
example, if lack of this mass term in the ghost particle propagator, the W-T identities shown in Eqs. (4.13), (5.14) 
and (5.15), which are derived from Eqs. (4.6) and (5.9) respectively, will have different expressions. In the lowest 
order approximation, these expressions can not be converted to the results which coincide with the Feynman rules. 
At this point, we can say, some previous works'^^l are not correct because these theories did not give the Lagrangian 
a ghost particle mass term in the general gauges. 

At present, the massless QCD has widely been recognized to be the candidate of the strong interaction theory 
and has been proved to be compatible with the present experiments. However, we think, the QCD with massive 
gluons would be more favorable to explain the strong interaction phenomenon, particularly, at the low energy domain 
because the massive gluon would make the force range more shorter than that caused by the massless gluon. As 
for the high energy and large momentum transfer phenomena, as seen from the massive gluon propagator, the gluon 
mass gives little influence on the theoretical result in this case so that the massive QCD could not conflict with the 
well-established results gained from the massless QCD in the high energy region. 
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VIII. APPENDIX: VERIFICATION OF THE SOLUTION TO THE W-T IDENTITY SATISFIED BY THE 

COUNTERTERMS 

In this appendix, it is shown that the solution given in Eq. (2.27) indeed satisfies the W-T identity written in Eq. 

(2.25) (or in Eq. (2.21)). Since the Ha in the first term of Eq. (2.27) are the gauge-invariant functional of the field 
functions A^, V' and if), according to the definition in Eq. (2.19), it is seen that^^^l 



SSa 5H^ I SSo SH^ _ SSo SH^ 



Scpi Sui Sui Scfii 5ui Sipi ( 

where the fact that the Ha are independent of the source variables and gauge-invariant, SHa ~ 0, has been noted. 
For the second term in Eq. (2.27), it is only needed to prove that the operator p{So) is nilpotent'^°l. Let us calculate 
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fjyOQ) r — Sipj ' Suj V Slp^ ' Su, ^ Sui ' Sipi I 
I SSq , _§_(SSa . I 5So . SF„ \ 

_ SSo_ . ( S^Sq . SFr^ I iSo . S^F„ 

I 5^ So , 5F^ _ SSo . S^Fn \ 

SujSui Sipi Sui SujSipi^ 
I SSo . ( S'\S„ , 5I;\ I dSo S'^ F„ /^2') 

I S^So , I 5So . S^F„ \ 
SipjSui Sipi Sui SipjS^pi' 
_ _S_fSSo_ . SSo_\ . SF\_ _|_ _5_ciSo . SSo_\ . jF„ 
(5(pi V Spj Suj y Sui Sui V Sip^ Suj / (Stpi 

I i5So . 5So . S'^F„ _ SSo . SSo . (5^F„ 
(5wj (5(pi SipjSui Sipj Sui SujS^i 

I SSq . dSo . i5^F„ I ilo . 5So . iS^Fn 
(5<^j <5(^i SujSui 5uj Sui SipjSipi' 

In the above, we assume that the field variables are commuting and the source variables are anticommuting for 

convenience of statement. With this assumption, It is easy to see that the third and fourth terms are cancelled with 
each other and the fifth and sixth terms themselves are identical to zero because |^ • 4^ and are symmetric 

tensors, while and |^ • are antisymmetric ones. Therefore, we have 

where {§0 * S) = has been noted. Similarly, it can be proved 

piS„-ifFn = 0. (A4) 
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X. FIGURE CAPTIONS 

Fig. (1): The one- loop gluon self-energy. The solid, wavy and dashed lines represent the free quark, gluon and 
ghost particle propagators respectively. 

Fig. (2): The one-loop ghost particle self-energy. The lines represent the same as in Fig. (1). 
Fig. (3): The one-loop ghost-gluon vertices . The lines mark the same as in Fig. (1). 

Fig. (4): The one- loop effective coupling constants aii(A) given by the timelike momentum space subtraction. The 
solid and dashed lines represent the coupling constants given by the massive QCD and the massless QCD, respectively. 
The dotted line denotes the coupling constant given in the minimal subtraction scheme. 

Fig. (5): The one-loop effective coupling constants aii{X) given by the spacelike momentum space subtraction. 
The solid and dashed lines represent the coupling constants given by the massive QCD and the massless QCD, 
respectively. The dotted line denotes the coupling constant given in the minimal subtraction scheme. The dashed- 
dotted line represents the one- loop effective coupling constants q;_r(A) given by the spacelike momentum subtraction 
for which the gluon mass is taken to be a smaller value. 
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